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ABSTRACT 


Pulse  Doppler  radar  systems  are  increasingly  being  used  for  airborne  maritime 
surveillance  roles,  but  current  models  for  sea  clutter  do  not  incorporate  the 
spectral  information  required  to  analyse  the  detection  performance  of  these 
systems  in  the  low  grazing  angle  maritime  environment.  This  report  uses  stan¬ 
dard  radar  sea  clutter  models,  and  spectral  information  from  wavetank  and 
open  ocean  measurements  to  deduce  the  amplitudes,  Doppler  frequencies,  and 
bandwidths  of  the  Bragg,  whitecap  and  spike  components,  which  together  com¬ 
prise  the  sea  clutter  return.  This  version  of  the  report  incorporates  corrections 
listed  in  the  Corrigenda. 
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A  Model  of  Low  Grazing  Angle  Sea  Clutter  for  Coherent 
Radar  Performance  Analysis 

Executive  Summary 

Airborne  pulse  Doppler  radar  systems  developed  for  air-to-air  applications  are  increasingly 
being  used  in  the  maritime  surveillance  role.  Because  of  the  relatively  low  peak  powers 
of  these  radars  compared  with  dedicated  maritime  surveillance  radars,  they  tend  to  use 
coherent  processing  to  enhance  the  signal-to-clutter-plus-noise  ratio  sufficiently  for  target 
detection.  In  analysing  the  performance  of  these  radars,  conventional  low  grazing  angle 
models  for  the  sea  clutter  are  inadequate  as  they  provide  only  the  mean  radar  cross  section 
for  horizontal  and  vertical  polarisation.  To  correctly  estimate  coherent  radar  performance, 
the  Doppler  characteristics  of  sea  clutter  are  required,  but  there  are  no  models  describing 
this  aspect  of  sea  clutter  in  the  open  literature.  This  work  seeks  to  develop  a  suitable 
Doppler  model  from  conventional  clutter  models  based  on  information  on  the  Doppler 
characteristics  of  sea  clutter  from  wavetank  and  ocean  data  collection. 

Three  conventional  clutter  models,  namely  the  GIT,  Hybrid  and  NRL  models,  which 
provide  estimates  of  the  mean  radar  cross  section  for  horizontal  and  vertical  polarisation 
as  a  function  of  radar  frequency,  sea  state,  grazing  angle,  and  wind  aspect  angle,  are 
documented.  They  are  then  modified  to  comply  with  recently  published  measurements 
of  the  variation  of  vertical  polarisation  clutter  with  aspect  angle  which  differs  from  the 
original  data  used  to  develop  these  models.  We  then  review  research  on  the  spectral 
characteristics  of  waves,  which  shows  that  sea  clutter  can  be  described  by  three  separate 
components,  Bragg  scattering,  whitecap  scattering  and  spike  scattering,  each  with  a  char¬ 
acteristic  Doppler  peak  frequency  and  spectral  bandwidth.  Then,  based  on  assumptions 
of  the  relative  amplitudes  of  these  components  as  a  function  of  wind  speed,  and  the  peak 
frequencies  and  bandwidths  obtained  from  wavetank  data,  we  describe  two  methods  for 
deducing  the  amplitudes  of  the  spectral  components  from  the  conventional  clutter  models. 

Charts  are  provided  of  the  Bragg,  whitecap  and  spike  terms  for  both  horizontal  and 
vertical  polarisations  as  a  function  of  wind  aspect  angle  for  sea  states  1  to  5,  determined 
by  the  two  methods,  for  each  of  the  conventional  clutter  models.  There  is  considerable 
similarity  of  form  between  the  clutter  models,  the  differences  mainly  resulting  from  the 
different  clutter  levels  predicted  by  the  models.  Which  clutter  model  is  selected  may 
depend  upon  the  requirements  of  the  study  for  which  the  spectrum  is  required.  The  two 
methods  have  slightly  different  variations  as  a  function  of  wind  aspect,  but  it  is  not  possible 
from  current  sea  clutter  data  to  select  one  as  being  preferred.  There  is  little  difference  in 
the  overall  spectra  from  the  three  clutter  models,  and  only  that  obtained  from  the  NRL 
model  are  presented  in  this  report.  A  comparison  of  the  results  with  the  limited  published 
open  ocean  spectra  data  indicates  that  the  model  described  will  be  a  useful  tool. 
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1.  page  5,  equation  13,  second  line,  changed  exponent  of  Hav  +  0.015  from  -0.24  to  -0.4. 

2.  page  5,  equation  13,  both  lines,  changed  if)  +  0.001  to  V’  +  0.0001. 
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1  Introduction 

Sea  clutter  has  been  studied  extensively  since  the  development  of  airborne  maritime 
surveillance  and  naval  radars  in  World  War  2.  In  the  version  of  the  radar  equation  for 
estimating  performance  of  these  radars  [7] ,  the  effect  of  sea  clutter  is  included  as  a  simple 
term  do,  the  radar  cross  section  per  unit  area  of  the  ocean  surface,  which  is  then  multiplied 
by  the  effective  area  of  the  radar  resolution  cell  to  provide  the  total  clutter  signal  against 
which  the  target  signal  has  to  be  detected.  What  is  not  apparent  in  the  radar  equation  is 
that  this  term  is  dependent  on  a  large  number  of  physical  parameters  of  the  ocean  envi¬ 
ronment,  the  wind  strength,  the  duration  of  the  wind,  and  the  fetch  (usually  these  three 
parameters  are  combined  as  a  sea  state  term),  the  wind  direction  relative  to  the  radar 
look  direction,  the  grazing  angle  of  the  radar  beam  with  respect  to  the  surface,  and  the 
presence  of  swell.  There  are  also  radar  system  parameters  which  can  impact  the  observed 
value  of  do,  the  radar  frequency,  the  pulse  repetition  frequency,  the  antenna  polarisation, 
and  the  resolution  cell  size.  Propagation  conditions  in  the  atmosphere,  especially  ducting, 
also  influence  the  observed  clutter  cross  section.  Nathanson  in  his  second  edition  of  Radar 
Design  Principles  [7]  cites  no  fewer  than  60  separate  clutter  measurement  studies  in  devel¬ 
oping  a  set  of  tables  of  the  mean  value  of  do  as  a  function  of  sea  state,  carrier  frequency, 
polarisation  and  grazing  angle.  Usually  the  radar  analyst  requires  clutter  cross  section 
data  under  conditions  differing  from  that  of  the  measured  data,  and  consequently  a  num¬ 
ber  of  formulas  have  been  developed  from  Nathanson’s  and  other  data  that  approximate 
the  functional  form  for  the  average  radar  cross  section  of  the  sea.  Notable  among  these  are 
the  Georgia  Institute  of  Technology  model  (GIT)  [6],  the  Technology  Service  Corporation 
model  (TSC)  [25],  and  a  Hybrid  model  [21],  The  US  Naval  Research  Laboratory  [5]  and 
the  UK  DERA  [3]  have  also  developed  models  for  their  own  applications. 

For  a  number  of  reasons,  coherent  pulse  Doppler  radars  developed  for  air-to-air  ap¬ 
plications  are  increasingly  being  used  for  maritime  surveillance  as  a  secondary  role,  and 
these  radars  tend  to  have  lower  peak  powers  than  dedicated  airborne  maritime  surveillance 
radars.  In  current  designs,  they  also  tend  to  have  longer  processed  pulse  lengths,  measured 
in  tenths  of  microseconds  rather  than  tens  of  nanoseconds.  Recourse  to  K  distribution1 
analysis  (see  for  example  Ward,  Baker  and  Watts  [28],  [1],  [29])  or  similar  is  generally  not 
required  as  conventional  Gaussian  statistics  are  adequate  to  describe  the  clutter  amplitude. 
These  radars  rely  upon  coherent  processing  to  improve  the  signal-to-clutter-plus-noise  ra¬ 
tio  sufficiently  for  target  detection,  and  thus  information  on  the  Doppler  characteristics 
of  the  sea  clutter  is  required  for  system  design  and  performance  estimation.  Models  to 
describe  the  effects  of  low  grazing  angle  sea  clutter  on  pulse  Doppler  radar  system  perfor¬ 
mance  over  the  range  of  parameters  of  operational  interest  currently  do  not  exist  in  the 
open  literature. 

Observations  of  sea  clutter  phenomena,  from  wavetank  experiments  and  direct  mea¬ 
surement  of  the  ocean  returns,  conclude  that  there  are  three  separate  scattering  mech¬ 
anisms  that  need  to  be  included  in  the  description  of  the  Doppler  characteristics  of  sea 
clutter.  These  mechanisms  are  resonant  Bragg  scattering  by  capillary  waves  which  are 
modulated  by  the  longer  gravity  wave  structure,  scattering  by  very  rough  whitecaps  pro¬ 
duced  by  breaking  waves,  and  specular  reflection  from  the  crest  of  the  wave  just  before 

XK  is  the  modified  Bessel  function  of  the  first  kind 
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it  breaks,  producing  sea  spikes  [26],  [8].  Both  amplitude  and  Doppler  characteristics  of 
the  scattering  mechanisms  are  required.  The  relative  significance  of  these  phenomena  de¬ 
pends  upon  the  incident  polarisation.  For  vertical  polarisation,  the  Bragg  scattering  term 
is  dominant,  whereas  with  horizontal  polarisation  this  term  is  only  weakly  backscattered. 
Scattering  by  whitecaps  appears  to  be  essentially  polarisation  independent,  whereas  the 
sea  spike  mechanism  produces  a  much  stronger  horizontal  polarisation  response  and  is  vir¬ 
tually  absent  with  vertical  polarisation.  There  is  controversy  in  the  literature  over  whether 
Bragg  scattering  occurs  at  microwave  frequencies,  see  for  example  articles  by  Skolnik  [24] 
and  Wetzel  [31].  The  US  Naval  Research  Laboratory  has  observed  that  at  very  high  reso¬ 
lution,  all  clutter  appears  spiky  even  at  low  sea  states,  and  that  there  is  not  the  continuum 
of  echo  signals  suggested  of  Bragg  scattering  models.  However  the  majority  of  papers  on 
spectral  characteristics  of  sea  clutter  refer  to  Bragg  scattering  as  being  the  contributor  of 
one  component  of  the  clutter,  often  the  most  significant,  and  we  shall  continue  with  this 
terminology  as  the  origin  of  this  signal  is  not  important  in  this  study. 

In  this  report  we  develop  a  technique  for  determining  the  amplitude  and  Doppler  char¬ 
acteristics  of  the  differing  components,  in  which  the  spectral  parameters  of  peak  frequency 
and  bandwidth  are  obtained  from  wavetank  measurements  of  the  Doppler  characteristics 
of  waves,  and  then  the  relative  amplitudes  of  the  spectral  components  are  deduced  from 
standard  models  for  no-  Comparison  are  made  of  the  predictions  obtained  with  the  GIT, 
Hybrid  and  NRL  clutter  models.  We  find  that  the  spectrum  of  open  ocean  waves  is  wider 
than  that  observed  in  wavetanks,  and  the  wavetank  models  need  to  be  multiplied  by  a 
scale  factor  to  agree  with  measurements  of  ocean  waves.  Although  the  emphasis  in  this 
report  is  on  the  development  of  the  clutter  spectral  model,  it  is  only  a  means  to  an  end; 
to  determine  the  probability  of  detection  of  a  target  in  sea  clutter  when  using  coherent 
processing,  knowledge  is  required  of  the  clutter  covariance  matrix,  with  which  the  im¬ 
provement  in  signal-to-clutter-plus-noise  ratio  as  a  result  of  the  coherent  processing  can 
be  determined.  Farina  et.  al.  [4]  have  described  a  technique  for  determining  this  improve¬ 
ment  factor  for  a  clutter  model  developed  by  Lee  et.  al.  [10],  and  it  is  simple  to  adapt 
to  the  spectral  clutter  model  that  is  developed  in  this  report.  On  the  supposition  that 
the  Doppler  components  are  statistically  independent,  it  is  relatively  straightforward  to 
calculate  the  clutter  improvement  factor  using  this  method. 

Plant  and  others  have  published  a  number  of  papers  in  the  literature  to  indicate  that 
a  theoretical  investigation  of  the  spectra  for  vertical  polarisation  sea  clutter  could  be  a 
profitable  avenue  of  investigation,  but  such  an  investigation  was  not  considered  likely  to 
be  completed  in  time  for  the  particular  requirement  for  which  this  clutter  model  was 
required.  Plant  et  al.  [19]  extended  a  composite  scattering  model  down  to  grazing  angles 
of  1°  by  incorporating  shadowing  by  the  waves  into  the  clutter  calculations,  and  showed 
that  (Jq(VV)  is  explained  by  the  model  over  the  whole  range  of  grazing  angles  to  zenith, 
but  measured  (Jq{HH )  tends  to  be  greater  than  the  model  predicts  for  grazing  angles 
less  than  about  45°.  In  an  earlier  paper,  Plant[15]  discussed  the  problem  of  modelling  of 
Doppler  effects  in  a  one  dimensional  wave,  and  showed  that  it  is  possible  to  qualitatively 
predict  the  observed  effects  of  the  Doppler  spectra  of  <jq  {HH)  and  no  (VV).  However  he 
did  not  pursue  this  matter  further  in  his  subsequent  papers. 

In  the  literature  on  sea  scatter,  the  angle  of  the  radar  beam  to  the  ocean  surface  is 
described  by  either  an  incidence  angle  or  a  grazing  angle,  the  former  being  more  commonly 
used  in  the  remote  sensing  literature.  The  incidence  angle  refers  to  the  angle  the  radar 
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beam  propagation  direction  makes  with  the  normal  to  the  ocean  surface,  whereas  the 
grazing  angle  is  the  angle  of  the  beam  to  the  (tangential)  plane  of  the  ocean  surface.  We 
will  use  the  latter  in  this  report. 

In  the  documents  referenced  in  this  report  there  is  no  uniformity  in  the  use  of  symbols 
for  the  input  parameters  to  the  models,  the  sea  state,  the  wind  speed  and  various  wave 
parameters,  the  grazing  angle,  and  the  wind  aspect,  so  where  necessary  the  notation  has 
been  changed  from  the  original  documentation  to  ensure  consistency  within  this  report.  In 
addition,  the  values  of  numerical  constants  have  been  changed  from  the  original  references 
so  that  across  all  the  models,  within  formulas  distances  are  expressed  in  metres,  velocities 
in  metres/second,  and  angles  in  radians. 


2  Amplitude  Characteristics  of  Sea  Clutter 

Most  of  the  early  work  on  sea  clutter  concentrates  on  the  amplitude  characteristics  of 
the  clutter,  and  several  differing  models  to  describe  the  average  radar  cross  section  of  sea 
clutter  have  been  developed.  The  ones  that  are  included  in  this  study  are  the  GIT,  Hybrid 
and  NRL  models.  Our  interest  is  in  fully  developed  seas,  primarily  sea  states  1  to  5,  in 
which  it  is  assumed  that  the  appropriate  wind  strength  has  been  blowing  for  the  required 
duration,  and  that  the  fetch  is  sufficient  for  the  seas  to  fully  develop.  The  sea  state, 
specifically  the  Douglas  sea  state,  is  an  integer  between  0  and  9,  each  of  which  describes 
a  range  of  sea  (and  by  implication  wind)  conditions.  Within  the  clutter  models  though, 
sea  state  is  treated  as  a  continuous  parameter,  from  which  a  wind  velocity  is  determined 
via  the  formula 

14  =  3.16  SS0-8  (1) 

Table  1  compares  the  wind  range  for  each  sea  state  with  the  wind  velocity  used  in  the 
models,  and  it  can  be  seen  that  the  models  calculate  the  clutter  in  the  upper  portion  of 
the  range  described  by  each  sea  state  index. 

Table  1:  Comparison  of  Douglas  sea  state  wind  ranges  with  sea  state  algorithm 


Sea  state 

wave  height  (m) 

wind  range  (m/s) 

model  wind  (m/s) 

1 

0-0.1 

0-3.6 

3.2 

2 

0. 1-0.5 

3. 6-6. 2 

5.5 

3 

0.5-1.25 

6. 2-8. 2 

7.6 

4 

1.25-2.5 

8. 2-9. 8 

9.6 

5 

2. 5-4.0 

9.8-11.8 

11.5 

For  a  fully  developed  sea,  there  is  a  simple  relationship  between  the  wind  speed  and  the 
parameters  which  describe  the  sea  surface.  In  the  GIT  model  the  following  relationships 
are  used 


Lav  =  0.08  SS2, 


(2) 


and 

14  =  8.67  hlA 


(3) 
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where  hav  is  the  average  wave  height  in  metres,  SS  is  the  sea  state,  and  Vw  is  the  velocity 
of  the  wind  in  metres/second.  The  Hybrid  model  also  uses 

ft-RMS  =  0.031  SS2  (4) 

where  /irms  is  the  RMS  wave  height.  In  describing  sea  conditions,  the  significant  wave 
height,  hi/ 3,  the  average  height  of  the  one  third  highest  waves,  is  the  primary  indicator  of 
surface  roughness,  and  the  NRL  model  uses  the  relationship 

hi/3  =  0.049  SS2'6  (5) 

though  it  does  not  appear  directly  in  their  equations.  Pidgeon  [14]  lists  a  number  of  other 
relationships  of  the  “average”  wave  parameters  for  a  given  wind  velocity  which  are  useful 
but  are  not  of  significance  for  the  research  in  this  report. 


2.1  The  GIT  Sea  Clutter  Model 


The  Georgia  Institute  of  Technology  empirical  sea  clutter  model  [6]  has  widespread  accep¬ 
tance  as  an  appropriate  model  for  the  mean  radar  cross  section  at  microwave  frequencies 
at  low  grazing  angle  when  there  is  no  anomalous  propagation.  It  originally  contained 
two  separate  expansions  covering  the  frequency  ranges  1  to  10  GHz,  and  10  to  100  GHz, 
but  subsequently  the  break  point  was  changed  to  12  GHz  to  include  all  of  X  band  in  the 
former.  The  basis  of  the  model  is  that  at  low  grazing  angles,  theoretical  considerations, 
supported  by  observations  obtained  in  non-ducting  conditions,  indicate  that  uq  varies  as 
the  fourth  power  of  the  grazing  angle,  ip,  and  above  a  critical  angle  flattens  off  in  a  plateau 
region,  before  increasing  rapidly  with  angles  near  zenith.  Observations  culled  to  retain 
only  non-ducting  conditions,  were  then  used  to  evaluate  the  parameters  in  the  model.  We 
are  concerned  with  only  the  lower  frequency  range  expansion  in  this  study. 

The  GIT  model  comprises  three  separate  terms,  a  sea  state  term,  Aw,  a  multipath 
term,  At,  and  a  sea  direction  term,  Au,  each  of  which  are  functions  of  several  variables. 
Data  used  in  the  development  of  the  model  indicated  that  the  sea  state  term  is  more  a 
function  of  wind  speed  Vw,  than  wave  height.  Plotting  ao  against  the  logarithm  of  the 
wind  speed  provided  approximately  straight  lines  of  the  form  (A  +  0.015) _0'4  for  low  wind 
speeds,  where  A  is  the  radar  wavelength  (in  metres).  Taking  account  of  the  saturation 
effect  at  high  speeds,  the  wind  speed  factor  became 


A 


W 


1-94  Vw 

1  +  Uo/15.4 


1.1  (A+0.015)-04 


(6) 


For  a  fully  developed  sea,  wind  speed  and  the  average  height  of  the  sea  are  interrelated  by 


Vw  =  8.67  hJA 


(7) 


where  Vw  is  in  meters  per  second  and  hav  is  in  metres. 

The  multipath  term  is  based  on  interference  between  a  direct  ray  and  a  (ocean)  surface 
reflected  ray,  the  amplitude  of  which  is  modified  by  the  surface  roughness.  A  surface 
reflection  coefficient  of  -1  was  assumed  for  all  grazing  angles,  though  this  is  strictly  only 
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accurate  for  horizontal  polarisation.  Scattering  from  a  Gaussian  surface  suggests  that  the 
roughness  coefficient  should  be 

J -  (8) 

where  ip  is  the  grazing  angle,  and  /iav  is  the  average  wave  height.  However,  Horst  et  al. 
observed  from  experimental  data  that  the  roughness  parameter  is  wavelength  dependent 
and  they  modified  Equation  8  to  read 


(14.4A  +  5.5)  ip  /iav 


(9) 


It  is  assumed  that  ip  is  sufficiently  small  that  sin  ip 
becomes 


Ai 


oy.4 

l  +  oy4 


ip.  The  interference  term  then 


(10) 


The  sea  direction  term  is  based  on  upwind-downwind  ratio  data,  and  is  dependent  on 
the  grazing  angle  ,  increasing  for  decreasing  grazing  angle.  The  dependence  on  wind  aspect 
<p  was  based  on  the  observation  that  <7o  for  HH  polarisation  in  the  crosswind  direction  is 
about  midway  in  decibel  value  between  the  upwind  and  downwind  values.  As  there  was 
no  data  to  define  the  functional  dependence,  a  sinusoidal  function  was  chosen, 

Au  =  exp[ 0.2  cos  <p(l  —  2.8 ip)(X  +  0.015)-0'4]  (11) 


for  both  polarisations,  where  <p  is  the  angle  between  the  radar  look  direction  and  the  wind 
direction,  having  the  value  0°  when  the  radar  is  looking  upwind. 

The  average  cross  section  per  unit  area  of  the  ocean  surface  for  horizontal  polarisation 
transmission  and  reception,  <JqH  ,  is  then 

aU11  =  3.9  x  10~6  A  ip0AAi  Au  Aw  (12) 

in  which  the  constant  term  has  been  experimentally  determined,  A  is  a  residual  wavelength 
dependency,  and  ip0A  accounts  for  the  effects  of  shadowing  at  low  angles.  The  vertical 
polarisation  cross  section,  <7q  v ,  differs  from  (JqH  by  a  multiplicative  factor,  the  form  of 
which  depends  whether  the  radar  frequency,  /RF,  is  greater  or  less  than  3  GHz.2 

f  9.33  (L/av  +  0.015)-0'24  A0'25  (ip  +  O.OOOl)0'29^,?^  /RF  >  3  GHz 
4V  =  {  (13) 

!  166.0  (iLav  +  0.015)-0’4  A0,87  (ip  +  O.OOOl)0’57^^  /RF  <  3  GHz 


2.2  Hybrid  Sea  Clutter  Model 

The  Hybrid  sea  clutter  model  was  first  described  in  a  limited  release  NATO  AAW  clut¬ 
ter  working  group  document  prepared  by  J.P.Reilly  [20]  of  Applied  Physics  Laboratory, 
Johns  Hopkins  University.  It  was  intended  to  provide  sea  clutter  estimates  under  specific 

2The  formulas  for  in  [6]  are  written  in  the  form  of  two  expressions  for  decibels  to  be  subtracted 
from  (To  in  dB.  These  are  converted  in  the  above  equations  to  multiplicative  terms  in  linear  units  for 
implementation  in  computer  software. 
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evaporation  duct  propagation  conditions,  but  at  the  time  of  release  the  propagation  study 
was  incomplete.  The  model  was  based  on  published  empirical  fits  to  clutter  data  that  in¬ 
cluded  some  ducting,  including  data  from  Barton  [2]  and  Nathanson  [7],  and  incorporated 
features  of  the  GIT  model.  The  frequency  range  of  the  data  base  used  was  from  0.5  to 
35  GHz.  Subsequently  Reilly  and  Dockery  [21]  showed  that  a  significant  source  of  vari¬ 
ability  in  sea  clutter  measurements  at  low  grazing  angles  can  be  attributed  to  variations  in 
the  atmospheric  refractive  index  changing  propagation  conditions  near  the  ocean  surface. 
Accounting  for  duct  induced  changes  in  grazing  angle,  and  in  space  propagation  factor, 
the  output  of  the  GIT  model  was  significantly  increased  and  produced  results  closer  to 
the  Hybrid  model. 

The  Hybrid  model  is  expressed  in  dB  and  has  the  form 

a0dB  =  do  (ref)  +  Kg  +  Ks  +  Kp  +  Kd  (14) 

in  which  cro(ref)  is  a  reference  reflectivity  applying  to  sea  state  5,  grazing  angle  0.1°, 
look  direction  upwind,  and  vertical  polarisation.  Kgi  Ks,  Kp,  and  Kd  are  adjustment 
factors,  to  accommodate  changes  from  the  reference  reflectivity  in  grazing  angle,  sea  state, 
polarisation  and  wind  direction.  The  reference  reflectivity  is 


22.4  log/RF  -  266.8  fRF  <  12.5  GHz 
3.25  log/RF  -  71.25  /RF  >  12.5  GHz 


(15) 


where  /rf  is  the  radar  frequency  in  Hertz. 
The  adjustment  factors  are 


1.  Kg  is  an  adjustment  for  changes  in  grazing  angle  'ip.  The  reference  grazing  angle  is 


Vv  =  0.1° 


(16) 


and  the  transition  angle3  is 


ipt  =  sin  1  (0.066  A/Lrms) 


where  Lrms  is  the  RMS  wave  height.  Then 
(a)  for  ipt  >  ipr 


Kg={ 


0 

20 

20 


log  Ip/lpr 

log  ipt/ ipr  T  10  log  ip/ipt 


for  ip  <  ipr 
for  ipr  <  ip  <  ipt 
for  ipt  <tp  <  30° 


(17) 


(18) 


(b)  for  ipt  <  ipr 


Kg  = 


0 

10  log  Ip /ipr 


for  Ip  <  1pr 
for  Ip  >  Ipr 


3In  [21],  the  numerical  factor  in  the  formula  is  0.66,  which  is  in  error. 


(19) 
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2.  The  sea  state  adjustment  used  by  Reilly  is 

Ks  =  5  (SS  -  5)  (20) 

which  places  the  lines  of  integer  sea-state  at  5  dB  intervals  at  ?/>  =  0.1°,  but  because 
the  transition  angle  increases  with  decreasing  sea-state,  the  curves  for  sea-states  1 
and  2  touch  for  horizontal  polarisation  and  cross  for  vertical  polarisation,  which  is 
clearly  not  correct.  To  counter  this  problem,  Ks  has  been  modified  in  this  work  to 
separate  out  the  curves  at  low  sea-state  as  follows4 

Ks  =  5  (SS  -  5)  +  (SS1~5)3  (21) 

with  minimal  effect  at  high  sea-state. 

3.  The  polarisation  adjustment,  Kp,  is  zero  for  vertical  polarisation.  For  horizontal 
polarisation  it  is  based  on  the  GIT  model5  and  has  the  value 

1.71n(/iav  +  0.015)  —  3.8  In  A 

-2.5  In  (?/>  + 0.0001)  -  22.2 
1.1  ln(/iav  +  0.015)  -  1.1A 

-  1.3  In  (V>  +  0.0001)  -  9.7 
1.41n(/iav)  —  3.4  In  A  —  1.3  Imj;  —  18.6 

4.  The  wind  direction  adjustment  Kd  is 

Kd=  ^2  + 1.7  log  (cos  <j)  1)  (23) 

where  is  the  angle  between  the  wind  direction  and  the  radar  look  direction,  set  to 
0°  when  the  radar  is  looking  upwind. 


/rf  <  3  GHz 

3  GHz  <  /RF  <  10  GHz 
/rf  >10  GHz 

(22) 


The  supplementary  equations  relating  wind  velocity,  RMS  wave  height  and  average 
wave  height  to  sea  state  used  in  this  model  are 


Vw  =  3.2  SS0'8 
h>RMS  =  0.031  SS2 
hav  =  0.08  SS2 


(24) 

(25) 

(26) 


2.3  The  NRL  Sea  Clutter  Model 

Gregers-Hansen  and  Mital  [5]  have  criticised  the  large  average  error  of  both  the  GIT 
and  Hybrid  models  when  compared  with  Nathanson’s  data.  In  particular  they  note  that 
the  GIT  model  considerably  underestimates  the  clutter  levels  at  low  sea-states.  They 

4In  [20]  Reilly  remarks  in  the  Commentary  and  Cautionary  Notes  that  the  approximation  of  5  dB  per 
sea  state  is  excessive  for  sea  states  above  6,  and  understated  for  sea  states  below  3. 

5  Note  that  all  the  numerical  factors  have  been  rounded  down  from  the  GIT  model 
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have  developed  an  empirical  model  seeking  to  minimise  the  absolute  average  error  in  dB 
from  Nathanson’s  data,  separately  optimising  the  parameters  for  horizontal  and  vertical 
polarisation,  of  the  form 

o-0dB  =  ci  +  c2log(sinV0  +  Q  +  C7(  1  +  SS)2+cs-+=9Ss  (27) 

1  +  c5  ip  +  C6  SS 

where  /rf  is  the  radar  frequency  in  GHz,  and  is  the  grazing  angle  in  radians6.  One  of 
the  claimed  advantages  of  this  model  is  that  it  includes  both  sea  states  0  and  6.  Note, 
however,  that  this  model  provides  only  the  upwind  radar  cross  section,  which  has  the 
greatest  reflectivity  of  all  wind  aspects.  The  parameters  of  Equation  27  are  listed  in 
Table  2. 

Table  2:  Parameters  for  the  NRL  clutter  model 


Parameter 

Polarisation 

Horizontal 

Vertical 

Cl 

-72.76 

-48.56 

C2 

21.11 

26.30 

C3 

24.78 

29.05 

c4 

281.7 

-29.70 

C5 

35.62 

60.56 

C6 

-0.02949 

0.04839 

C7 

26.19 

21.37 

C8 

5.354 

4.278 

C9 

0.05031 

0.04623 

Gregers-Hansen  and  Mital  provided  a  table  comparing  the  average  absolute  errors  of 
the  GIT,  Hybrid,  and  their  model  when  compared  with  Nathanson’s  data,  and  is  repro¬ 
duced  below 

Table  3:  RMS  errors  of  clutter  models  with  Nathanson’s  data 


Model 

Polarisation 

Horizontal 

Vertical 

GIT 

Hybrid 

NRL 

13.7  dB 

14.5  dB 

2.1  dB 

14.1  dB 

8.7  dB 
2.0  dB 

2.4  Revision  of  Angular  Dependence  of  Scattering 

In  the  GIT  sea  clutter  model,  the  aspect  dependence  is  based  on  upwind-downwind  ratio 
data  which  was  a  function  of  the  incidence  angle,  increasing  as  the  angle  decreases,  and 

6In  the  original  paper  ip  was  expressed  in  degrees,  but  in  the  parameter  table  the  constants  C4,  C5  and 
C8  have  been  multiplied  by  180/7T  so  that  it  can  be  expressed  in  radians. 
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proportional  to  A~0'4  at  ip  =  0.  For  HH  polarisation,  the  crosswind  RCS  was  observed  to  be 
about  midway  in  dB  between  the  upwind  and  downwind  RCS,  and  a  simple  cosine  function 
was  chosen  for  the  angular  dependence.  This  same  model  was  assumed  appropriate  for 
VV  polarisation.  Converting  Equation  11  to  dB  this  angular  dependence  becomes 

AudB  =  0.8686 cos </>(l  -  2.8^)(A+ 0.015)-0'4  (28) 

Note  that  there  is  no  dependence  on  the  sea  state  in  this  expression.  In  the  Hybrid  model, 
the  angular  dependence  in  dB  is  given  by  Equation  23,  and  is  independent  of  both  sea 
state  and  grazing  angle,  though  it  is  modified  by  the  wavelength. 

At  high  grazing  angles  both  HH  and  VV  polarisation  exhibit  a  different  angular  be¬ 
haviour.  The  RCS  exhibits  a  local  maximum  at  both  upwind  and  downwind  look  direc¬ 
tions,  upwind  being  the  greater,  with  minima  in  the  crosswind  directions,  i.e.  there  is 
a  second  harmonic  dependence  of  the  RCS  on  the  aspect  angle.  As  the  grazing  angle  is 
reduced,  there  is  point,  above  10°,  at  which  the  angular  behaviour  of  horizontal  polar¬ 
isation  changes  to  the  fundamental,  as  used  in  the  GIT  model.  However,  in  a  series  of 
measurements  in  the  South  China  Sea  in  2005,  and  off  the  New  Jersey  USA  coast  in  2006 
Plant  et  al.  [16]  found  that  at  low  grazing  angle,  vertical  polarisation  scatter  has  a  second 
harmonic  dependence  and  thus  is  at  variance  with  the  GIT  model. 

Plant  published  a  chart  of  vertical  polarisation  measurements  at  1°  grazing  angle  and 
a  wavelength  of  0.03  nr,  for  wind  speeds  of  4,  6,  and  8  m/s,  to  each  of  which  we  have  fitted 
a  three  term  cosine  series.  In  this  data,  the  variation  in  cross  section  with  aspect  increases 
with  increasing  wind  velocity.  However,  these  curves  have  been  averaged  to  provide  a 
result  independent  of  wind  velocity  to  comply  with  the  precepts  of  the  GIT  model;  the 
average  is 

Au  dB  =  1-67  cos0  +  3.15cos(2  (f>)  (29) 

If  the  dependence  of  the  GIT  model  on  t/>  and  A  is  retained,  but  we  otherwise  comply  with 
Plant’s  measurements,  then  a  revised  angular  dependence  model  for  VV  polarisation  takes 
the  form 


Al  =  exp  [(0.2206 cos(2</) +  0.1171  cos </>)  (1  -  2.8t/>)  (A +  0.015)-0'4]  (30) 

Plant  also  provides  measurements  of  the  Doppler  shifts  of  vertical  polarisation  corrected 
to  line  of  sight  velocities,  which  have  a  cosinusoidal  dependence7. 

Plant’s  measurements  of  horizontal  polarisation  at  a  grazing  angle  of  1.7°,  and  A  = 
0.03  nr  suggest  that  the  angular  dependence  should  be 

A^  =  4.43  cos  <f>  (31) 

Thus  the  revised  model  for  A l  in  the  GIT  model  becomes 

Al  =  exp  [0.3219 cos <£(1  -2.8^)  (A+  0.015)-0'4]  (32) 

In  the  Hybrid  model,  the  aspect  dependency  factor  is 

Kd  =  (0.3  —  1.71ogA)(cos0  —  1)  (33) 

7  Lamont-Smith  [9]  did  not  obtain  cosinusoidal  dependence  for  vertical  polarisation  in  his  measurements 
from  Portland  Bill,  only  for  horizontal  polarisation. 
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which  has  an  average  value  —(0.3  —  1.71ogA)  at  (j)  =  n/2.  If  we  scale  Plant’s  result  for 
vertical  polarisation  according  to  the  dependence  of  K d  on  the  radar  wavelength,  then  we 
have  the  result 


K%  =  (0.3-  1.71ogA)(0.578cos<£  +  1.09 cos(2 0)  -  1.0)  (34) 

The  result  for  horizontal  polarisation  is 

K%  =  (0.3  -  1.7 logA)(1.54cosc/>  -  1.0)  (35) 

Note  that  the  variation  in  amplitude  of  horizontal  polarisation  using  Plant’s  measurements 
is  somewhat  greater  than  in  the  original  Hybrid  model. 

The  NRL  model  does  not  include  the  effect  of  the  wind  aspect  on  the  radar  cross 
section,  and  does  not  include  guidance  on  how  the  other  parameters  affect  the  aspect 
dependency.  In  the  absence  of  any  other  information  we  shall  assume  a  functional  form 
identical  to  that  used  with  the  Hybrid  model. 

Figures  1  and  2  show  the  horizontal  and  vertical  polarisation  predictions  of  the  revised 
models  as  a  function  of  the  grazing  angle  for  sea  states  1  through  5  when  looking  upwind 
at  a  frequency  of  9300  MHZ.  The  GIT  and  NRL  models  are  comparable  at  high  sea  states 
and  large  grazing  angles,  but  at  low  grazing  angles  especially  for  sea  state  3  and  lower, 
the  predictions  of  the  GIT  model  fall  away  much  more  rapidly  than  those  of  the  NRL 
model.  The  Hybrid  model  gives  considerably  greater  predictions  than  the  other  two  at 
high  grazing  angles  for  all  sea  states.  In  its  original  form,  the  separation  between  curves 
at  low  sea  states  is  in  error  as  the  curves  for  sea  states  1  and  2  touch  for  horizontal 
polarisation  and  cross  in  vertical  polarisation. 

Figure  3  shows  the  horizontal  and  vertical  polarisation  predictions  as  a  function  of  wind 
aspect  at  a  grazing  angle  of  2°  for  the  GIT  model.  The  charts  for  the  angular  dependency 
of  the  Hybrid  and  NRL  models  are  similar,  the  variations  in  magnitude  are  the  same,  the 
only  difference  being  that  the  spacing  of  the  sea-state  lines  corresponds  to  the  different 
values  for  the  mean  radar  cross  section  predicted  by  these  models. 


3  Doppler  Characteristics  of  Sea  Clutter 

Notwithstanding  that  most  maritime  radars  in  the  1950s  and  1960s  used  incoherent  pro¬ 
cessing,  during  this  period  there  were  several  investigations  of  the  Doppler  characteristics 
of  sea  clutter.  The  most  significant  of  these  was  that  of  Pidgeon  [14],  in  which  he  took 
C  band  measurements  at  small  depression  angles  from  a  cliff  top  in  Puerto  Rico,  and 
subsequently  at  Cape  Cod,  Massachusetts.  He  observed  that  the  mean  Doppler  shift  for 
horizontal  polarisation  is  2  to  4  times  as  great  as  that  for  vertical  polarisation  in  the  same 
wind  and  wave  conditions.  He  concluded  that  the  vertical  polarisation  Doppler  shift  is 
directly  related  to  the  orbital  velocity  of  the  gravity  waves  whereas  the  horizontal  po¬ 
larisation  Doppler  shift  is  directly  related  to  the  motion  of  the  surface  layer  of  the  sea. 
For  both  polarisations,  the  Doppler  shift  was  directly  dependent  on  the  wave  height,  and 
varied  with  the  cosine  of  the  angle  between  the  wind  and  wave  direction  and  the  radar 
look  direction.  Also  he  reported  that  the  spectrum  bandwidth  was  relatively  independent 
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Grazing  angle  (degrees) 

(c)  modified  NRL  model 

Figure  1:  Predictions  of  a^H  (m2/m2),  of  the  modified  models,  expressed  in  dB,  as  a 
function  of  grazing  angle,  when  looking  upwind,  for  sea  states  1  to  5,  at  a  frequency  of 
9300  MHz 
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Grazing  angle  (degrees) 

(c)  modified  NRL  model 

Figure  2:  Predictions  of  a (m2/m2),  of  the  modified  models,  expressed  in  dB,  as  a 
function  of  grazing  angle,  when  looking  upwind,  for  sea  states  1  to  5,  at  a  frequency  of 
9300  MHz 
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Wind  aspect  (degrees) 

(a)  horizontal  polarisation 


Wind  aspect  (degrees) 

(b)  vertical  polarisation 

Figure  3:  Variation  with  wind  aspect  of  and  Oq  v,  expressed  in  dB,  for  the  modified 
GIT  model,  at  a  grazing  angle  of  2°  for  sea  states  1  to  5,  and  a  frequency  of  9300  MHz 


of  the  viewing  angle,  but  directly  dependent  on  the  mean  Doppler  shift  looking  into  the 
wind  and  waves. 

With  the  greater  use  of  Doppler  processing  in  most  modern  defence  radars  since  the 
mid  1980s,  there  has  been  considerably  greater  interest  in  the  Doppler  properties  of  sea 
clutter,  to  improve  the  detection  characteristics  of  targets  against  sea  clutter  with  this 
class  of  radars.  There  are  several  strands  of  research  that  we  will  follow,  notably  the  work 
described  in  the  papers  of  Lee  et  ah,  Walker,  Rozenberg,  and  Lamont-Smith. 


3.1  Basis  Function  Decomposition  of  Sea  Clutter  Spectra 

In  1991,  Lee  et  al,  [11]  participated  in  a  series  of  dual  polarisation  X  band  measurements 
of  ocean  spectra  at  Loch  Linnhe  and  the  Sound  of  Sleat  in  the  west  coast  of  Scotland, 
to  investigate  the  physics  of  microwave  backscattering  from  the  ocean  surface.  A  key  re- 
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suit  of  the  investigations  was  that  they  were  able  to  differentiate  Bragg  from  non-Bragg 
scattering,  and  to  resolve  the  difference  between  the  Doppler  structure  of  vertical  and 
horizontal  polarisations.  They  also  reported  instances  when  the  amplitude  of  the  horizon¬ 
tal  polarisation  return  exceeds  that  of  vertical  polarisation.  From  their  observations  of 
spiking  events,  they  concluded  that  spiking  events  can  mostly  be  related  to  wave  breaking 
events,  but  can  be  observed  at  other  times.  They  suggested  that  a  broader  definition  of 
spiking  events  was  required  to  describe  all  the  observed  phenomena  than  simply  stating 
that  it  originates  from  breaking  wave  events. 


Based  on  the  observation  that  differing  phenomena  accounted  for  the  spectral  com¬ 
ponents  of  sea  clutter,  Lee  et  al.  [10],  [12]  showed  that  the  clutter  spectra  could  be  de¬ 
composed  into  a  combination  of  Gaussian,  Lorentzian  and  Voigtian  basis  functions.  The 
Gaussian 


'M/) 


777 — ^exP 
WgVtt 


C f-fG ?' 
WG2  . 


(36) 


was  used  to  describe  Doppler  broadening  arising  from  a  spread  of  scatterer  speeds,  where 
fc  is  the  peak  frequency  of  the  spectrum,  and  Wq  is  the  spectral  width.  This  is  typi¬ 
cally  used  to  represent  Bragg  scattering.  The  Lorentzian  was  used  to  describe  lifetime  or 
damping-dominated  broadening  due  to  finite  scatterer  lifetime,  and  is  of  the  form 


*L(f) 


r/2vr2 

(/^)2+(02 


(37) 


In  this  expression,  fi  is  the  peak  Doppler  frequency  of  the  Lorentzian  spectrum,  and  T-1 
is  the  characteristic  scatterer  lifetime.  The  Voigtian  is  appropriate  if  the  scattering  arises 
from  both  lifetime  dominated  and  Gaussian  scattering  mechanisms,  and  is  a  convolution 
of  the  two  previous  profiles, 


H(a,  u) 


a 

IT 


exp(~  y2) 
(u  —  y )2  +  a2 


dy 


(38) 


in  which  a  =  (/  —  fv)/Wc,  fv  is  the  centre  of  the  Voigt  profile,  and  a  =  F/27tIUg  is  the 
ratio  of  the  Lorentzian  full  width  half  maximum  to  the  full  width  at  the  1/e  level  of  the 
Gaussian. 

Using  the  results  of  the  trials  at  Sound  of  Sleat  referred  to  above,  DRA  data  from  Loch 
Linnhe  that  was  collected  in  1989,  and  Pidgeon’s  C  band  data,  they  showed  that  VV  data 
can  be  represented  by  a  Gaussian  profile  with  a  peak  at  the  Bragg  resonant  frequency,  and 
a  Voigtian  contributing  to  the  higher  frequency  wing  of  the  power  spectral  density.  For 
HH  data,  the  best  representation  was  with  a  Lorentzian  for  the  principal  spectral  peak, 
and  a  Voigtian  contributing  to  both  the  upper  and  lower  wings  of  the  spectrum.  Some 
examples  of  the  fit  of  these  functions  to  the  experimental  data  from  [12]  are  reproduced 
in  Figure  4.  Note  that  in  these  charts  goodness  of  fit  of  the  functions  is  judged  on  a  linear 
scale. 

Lee  et  al.  [13]  subsequently  demonstrated  that  the  basis  functions  could  be  applied  to 
laboratory  data  (i.e.  wavetank  recordings)  as  well  as  ocean  data  to  describe  the  observed 
spectra. 
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60  80  100  120  HO  160  180  200  220  240 

frequency,  Hz 


(a)  TRW  data,  vertical  polarisation,  upwind,  wind  speed  9  m/s;  Gaussian,  dashed;  Voigtian, 
dotted;  sum  of  components,  solid 


60  8C  100  120  140  160  180  200  220  2  40 

frequency.  Hz 

(b)  TRW  data,  horizontal  polarisation,  upwind,  wind  speed  9m/s;  Lorentzian,  dashed;  Voigtian, 
dotted;  sum  of  components,  solid 


0  20  40  60  80  100  120  U0 

frequency.  Hz 


(c)  DRA  data,  vertical  polarisation,  upwind,  wind  speed  4  m/s;  Gaussian,  dashed;  Voigtian,  dotted; 
sum  of  components,  solid 

Figure  4:  Examples  of  fitting  Gaussian,  Lorentzian,  and  Voigtian  basis  functions  to 
recorded  sea  clutter  spectra  (from  Lee  et  al.  [12]) 
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3.2  Wavetank  Investigations  of  the  Doppler  Spectra  of  Waves 

At  about  the  same  time  period  as  the  work  of  Lee  et  ah,  Rozenberg,  Quigley  and  Melville 
[22],  [23]  measured  the  spectra  of  waves  at  14  GHz  (A  =  0.021  m)  in  a  wave  tank  operated 
by  the  Scripps  Institution  of  Oceanography  in  La  Jolla  California.  The  instrumentation 
enabled  simultaneous  measurement  of  the  HH  and  VV  polarisation  signals,  and  the  am¬ 
plitude  spectrum  of  the  waves.  Wind  speeds  were  in  the  range  2  to  12  m/s,  (roughly  4  to 
24  knots)  and  measurements  were  taken  both  upwind  and  downwind  at  grazing  angles  of 
6°,  8°,  and  12°.  The  main  results  of  interest  to  us  are  a  series  of  4  graphs  showing  the  peak 
of  the  Doppler  spectrum,  and  the  spectral  width  versus  wind  speed  for  both  polarisations 
as  a  function  of  wind  speed. 

Looking  downwind,  it  is  difficult  to  distinguish  the  peak  Doppler  frequencies  /  of 
HH  and  VV  polarisation,  and  there  is  negligible  variation  with  the  incidence  angle.  The 
observed  Dopplers  increase  approximately  linearly  with  the  wind  velocity.  Looking  upwind 
the  peak  Doppler  for  VV  polarisation  is  very  similar  to  downwind.  However,  the  Doppler 
for  HH  polarisation  is  considerably  greater  and  shows  a  slight  flattening  at  high  wind 
velocities.  Again  there  is  not  significant  variation  with  the  incidence  angle. 

The  downwind  bandwidths,  W,  of  vertical  and  horizontal  polarisation  are  also  difficult 
to  distinguish,  and  do  not  vary  significantly  with  grazing  angle.  Looking  upwind,  the 
vertical  polarisation  bandwidth  is  roughly  linear  with  frequency  up  to  the  limit  of  the 
wind  velocity  measurements,  but  the  horizontal  polarisation  bandwidth,  although  linearly 
proportional  for  small  wind  velocities,  plateaus  to  an  approximately  constant  value  at 
about  5  m/s. 

The  data  from  these  graphs  has  been  digitised  and  appropriate  functions  fitted  using 
the  FindFit  function  of  Mathematica.  The  data  on  the  peak  frequencies  of  the  spectrum, 
and  the  functions  fitted  to  the  data  are  shown  in  Figure  5,  and  on  the  bandwidth  are 
shown  in  Figure  6.  No  attempt  has  been  made  to  distinguish  the  data  points  from  grazing 
angles  of  6°,  8°,  and  12°  in  these  figures.  The  mathematical  approximations  displayed  in 
the  figures  are 

•  Upwind,  vertical  polarisation  /  =  17.36  +  10.59 U0'29  +  0.01531/3'05 

•  Upwind,  horizontal  polarisation  /  =  —39.43  +  57.48 VU  —  5.69 U 

•  Downwind,  both  polarisations  /  =  22.83  +  2.84 U 

•  Upwind  vertical  polarisation  W  =  5.28 U 

•  Upwind,  horizontal  polarisation  W  =  Min(6.15  U,  36.7) 

•  Downwind,  both  polarisations,  W  =  3.92 U 

where  U  is  the  wind  speed  in  metres/second,  and  Min  is  the  minimum  value  of  the  enclosed 
parameters.  Note  that  these  results  need  to  be  appropriately  scaled  to  provide  estimates 
of  the  position  of  the  spectral  peak  and  bandwidths  at  radar  frequencies  other  than  that 
used  in  these  measurements. 
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U  (m/s) 

(a)  vertical  polarisation,  upwind 


U  (m/s) 

(b)  horizontal  polarisation,  upwind 


U  (m/s) 

(c)  downwind,  both  polarisations 


Figure  5:  Peak  frequency  of  Doppler  spectrum  at  If  GHz  as  a  function  of  windspeed, 
data  points  from  [22],  solid  line  -  best  fit  of  mathematical  model 
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U  (m/s) 

(a)  vertical  polarisation,  upwind 


U  (m/s) 

(b)  horizontal  polarisation,  upwind 


Figure  6:  Bandwidth  of  Doppler  spectrum  at  If  GHz  as  a  function  of  windspeed,  data 
points  from  [22],  solid  line  -  best  fit  of  mathematical  model 
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3.3  Doppler  Measurements  of  Ocean  Waves 

The  papers  of  Pidgeon  and  of  Valenzuela  on  the  measurement  of  the  spectra  of  ocean  sea 
clutter,  while  instructive  on  the  characteristics  of  the  spectra,  do  not  contain  sufficient 
detail  to  develop  a  model.  However  in  2000,  Walker  [27]  published  a  three  component 
model  to  describe  wave  tank  measurements  similar  to  that  of  Lee  et  al.  above,  and  then 
applied  this  model  to  ocean  data  recorded  in  1996  from  a  clifftop  known  as  Portland  Bill 
on  the  south  coast  of  England  [26].  From  the  wavetank  work,  Walker  identified  three 
different  mechanisms  that  contributed  to  the  scattered  signal 

1.  Bragg  scattering,  arising  from  scattering  by  the  capillary  waves  riding  on  the  longer 
gravitational  waves,  both  polarisations  peaking  at  the  same  frequency,  and  being 
significantly  greater  in  magnitude  for  VV  polarisation  than  for  HH  polarisation 

2.  scattering  from  rough  whitecaps  of  broken  waves,  of  roughly  equal  magnitude  for 
both  polarisations,  and  speed  significantly  higher  than  the  Bragg  scattering 

3.  specular  scattering  from  the  crest  of  waves,  termed  spikes,  and  prominent  in  HH 
polarisation,  but  not  significant  in  VV,  with  a  Doppler  close  to  that  of  the  whitecaps 


Each  component  is  assigned  a  Gaussian  distribution.  For  vertical  polarisation,  the  model 
is  of  the  form 


,T,  ff\  r>v  f-(f-fB)2 

«V(/)  =  PB  «P  [  Wb2 

whereas  for  horizontal  polarisation 
*H(f)  =  Pb  exp  (  )  +Pw  exp  ( 


+  Pw  exp 

(/  ~  fw)2' 
Ww2  , 


-(/  ~  fw)2\ 

Ww2  ) 


+Ps  exp 


V  VF52  ) 


(39) 


(40) 


In  these  expressions,  subscripts  B,  w,  and  S  refer  respectively  to  the  Bragg,  whitecap  and 
spike  components.  Thus  /  is  the  radar  Doppler  frequency,  fs  and  fw  are  the  Doppler 
frequencies  corresponding  to  the  phase  velocity  of  the  Bragg  resonant  and  dominant  gravity 
waves  giving  rise  to  the  whitecaps,  and  Wb  and  Ww  are  the  widths  of  the  corresponding 
spectra.  Pj,  and  Pjf  are  the  amplitudes  of  the  vertical  and  horizontal  Bragg  components, 
whereas  Pw  and  Ps  are  the  whitecap  and  spike  components.  Note  that  the  whitecap 
component  is  the  same  in  both  polarisations,  and  that  the  spike  component  has  the  same 
Doppler  peak  as  the  whitecap  component,  though  a  different  spectral  width.  In  comparison 
with  Lee  et  al.  above,  the  use  of  Gaussian  functions  for  each  component  provides  very 
good  correlation  with  measurements  to  low  signal  levels  when  compared  on  a  logarithmic 
(dB)  scale. 

Two  sets  of  data  acquired  at  Portland  Bill  on  separate  days  are  presented  in  the  paper. 
The  radar  operated  at  a  centre  frequency  of  9.75  GHz,  at  64  m  above  mean  sea  level,  with 
a  grazing  angle  of  3.6°.  The  data  are  averaged  over  a  time  period  of  30  seconds.  The  first 
set  is  with  the  radar  looking  directly  into  the  wind,  with  a  windspeed  of  10  m/s,  gusting  to 
15  m/s.  The  second  set  is  with  the  radar  looking  downwind,  with  a  windspeed  at  7  m/s, 
gusting  to  10  m/s.  The  difficulty  with  the  measurements  is  that  detailed  information  on  the 
surface  could  not  be  obtained,  which  meant  that  many  surface  parameters  were  unknown 
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Figure  7;  Three  component  model  applied  to  upwind  data,  solid-  VV  model,  dashed-HH 
model  (from  Walker  [26]) 


Figure  8:  Three  component  model  applied  to  downwind  data,  solid-  VV  model,  dashed-HH 
model  (from  Walker  [26]) 


and  had  to  be  inferred  from  the  data.  In  the  upwind  VV  data,  the  Bragg  component  is 
clearly  visible,  and  the  whitecap  component  at  a  higher  frequency  is  only  slightly  lower 
in  amplitude.  In  the  HH  data,  the  Bragg  component  is  evident  as  a  slight  bump  in  the 
curve.  The  HH  peak  exceeds  the  VV  curve  at  this  point  indicating  a  contribution  from 
both  the  whitecap  and  spike  components  in  these  measurements.  The  downwind  spectra 
look  quite  different,  without  any  clear  peak  separation,  and  nowhere  does  HH  significantly 
exceed  VV.  The  Bragg  peaks  provide  the  greater  components  of  the  spectrum  with  the 
whitecap  component  much  smaller.  The  spike  term  does  not  contribute  to  the  downwind 
measurements.  Upwind  P\y  is  typically  4  to  10  dB  below  Pg ,  and  downwind  about 
14  dB  below  By  [26],  [27].  The  model  parameters  Walker  extracted  from  the  data  using 
a  Powell  minimisation  algorithm  are  listed  in  Table  4.  Figures  7  and  8  show  the  level 
of  agreement  with  the  data  of  the  three  component  model  for  the  upwind  and  downwind 
data  respectively  for  these  parameters. 
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Table  4 •'  Model  Parameters  for  Upwind  and  Downwind  Data  from  Portland  Bill  [26] 


WfflSk 1 

W3M3M 

MHz) 

M  Hz) 

WB  (Hz) 

Ww(  Hz) 

IBs  (Hz) 

Upwind 

-15.1  dB 

11.5  dB 

12.8  dB 

54 

178 

47.5 

66.0 

74.5 

Downwind 

-14.0  dB 

0.2  dB 

- 

-15.3 

-63.9 

40.4 

74.1 

- 

4  Constructing  a  Doppler  Model  of  Sea  Clutter 

There  is  no  published  model  which  satisfies  our  requirement  for  estimating  the  performance 
of  maritime  surveillance  modes  of  coherent  pulse  Doppler  radars,  but  there  are  many 
different  pieces  of  information  that  when  combined  will  allow  us  to  proceed  with  this  task 
in  a  qualified  way.  We  make  the  following  remarks 

1.  Walker’s  model  using  Gaussian  functions  appears  to  be  quite  appropriate  for  de¬ 
scribing  the  spectral  characteristics  of  the  ocean  conditions  of  interest;  it  provides 
good  results  over  a  significant  dynamic  range,  and  is  simpler  to  implement  than  the 
model  of  Lee  et  al.  A  key  feature  of  this  model  is  that  the  whitecap  component  is 
the  same  magnitude  for  both  polarisations. 

2.  Rozenburg  did  not  distinguish  the  origins  of  the  different  spectral  components  in  his 
wave  tank  measurements,  but  we  might  reasonably  conclude  from  applying  Walker’s 
conclusions  to  his  charts  that 

(a)  in  the  upwind  vertical  polarisation  measurements,  the  peak  in  the  wave  tank 
spectra  will  be  due  to  the  Bragg  component,  and  the  width  essentially  the  Bragg 
width 

(b)  in  the  upwind  horizontal  polarisation  measurements,  the  peak  in  the  wave  tank 
spectra  will  be  essentially  due  to  the  whitecap  component,  and  the  width  like¬ 
wise 

(c)  in  the  downwind  vertical  polarisation  measurements,  the  peak  in  the  wave  tank 
spectra  will  be  due  to  the  Bragg  component,  and  the  width  essentially  the  Bragg 
width 

(d)  in  the  downwind  horizontal  polarisation  measurements,  the  peak  in  the  wave- 
tank  spectra  will  be  due  to  the  Bragg  component,  but  the  width  will  be  broad¬ 
ened  due  to  the  presence  of  the  almost  equal  amplitude  whitecap  component 

3.  We  now  need  to  establish  the  positions  and  beamwidths  of  the  Bragg  and  whitecap 
terms  as  a  function  of  aspect  angle. 

(a)  Scaling  Rozenburg’s  measurements  of  the  Bragg  spectrum  peak  for  the  radar 
wavelength  A  (m)  of  interest,  from  the  upwind  measurements  we  have 

/g  =  (17.36  +  10.59  M29  +  0.0  1  53  M05)  0.021/A  (41) 

whereas  from  the  downwind  measurements 

/|  =  (22.83  +  2.84  U)  0.021  /A  (42) 
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where  superscripts  u.  and  d  indicate  upwind  and  downwind  directions.  At 
broadside  the  Doppler  of  the  Bragg  component  should  be  0,  so  an  appropriate 
function  for  the  angular  dependence  of  the  Doppler  might  be 


Ib 


fg  cos  4>,  0  <  <f>  <  7t/2, 

fgcosfi,  vr/2  <  4>  <  7r 


(43) 


Note  that  this  correctly  accounts  for  the  sign  of  the  Doppler  for  downwind 
cases. 

(b)  For  the  whitecap  term  from  the  upwind  measurements  the  peak  Doppler  is 

/&  =  (-39.43  +  57.48  VU  -  5.69 Uj  0.021/A  (44) 

and  from  the  downwind  measurements 


/((,  =  (22.83  +  2.84  U)  0.021/A 


(45) 


Again  zero  Doppler  is  expected  at  broadside,  so  an  appropriate  angular  depen¬ 
dence  is 

fw  cos  0,  0  <  (/>  <  7t/2, 


fw  = 


fw  cos  0;  n/2<(j)<TT 


(46) 


(c)  The  bandwidth  of  the  Bragg  component  follows  from  the  upwind  vertical  po¬ 
larisation  measurements 

Wl  =  5.28  U 


and  for  downwind 

W%  =  3.92  U 

Combining  these  results  with  a  sinusoidal  transition  provides 


WB  = 


w%  +  w£  iTh  -  wi 


+ 


2  2 
=  (4.6  +  0.68  U cos  <f>)  0.021/ A 


cos  <f> 


(47) 

(48) 

(49) 


(d)  The  bandwidth  of  the  whitecap  component  follows  from  the  horizontal  polari¬ 
sation  measurements.  Upwind  we  have 


Ww  = 


U  <  5.97 
U  >  5.97 


(50) 


and  downwind, 

W(v  =  3.92  U  (51) 

Combining  these  results  in  a  smooth  transition  with  angle 


Ww 


(5.035  +  1.115  cos 4>)  U 0.021/A  U  <  5.97 

(18.35  +  1.96  U  +  (18.35  -  1.96  U)  cos</>)  0.021/A  U  >  5.97 


(52) 
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(e)  Walker  records  that  the  spectral  width  of  the  spike  component  is  slightly  greater 
than  the  whitecap  component  by  a  factor  of  1.13  at  a  windspeed  of  10  m/s. 
Hence  we  shall  assume 

IDs  =  1-13  Ww  (53) 

(f)  We  note  at  this  stage  that  it  may  be  necessary  to  scale  the  above  frequency 
and  bandwidth  estimates  as  they  are  typically  much  greater  in  recorded  ocean 
data  than  in  data  obtained  from  wavetank  measurements. 

4.  We  have  been  discussing  spectral  power  densities,  but  clutter  models  such  as  the 
GIT  model  supply  average  cross  sections.  Integrating  Walker’s  spectral  model  over 
frequency  we  find  that 


(Tyy  =  s/tt  (ID b  Pb  +  Ww  Pw)  (54) 

and 

&hh  =  {Ww  Pb  +  Ww  Pw  +  Ws  Ps)  (55) 

Since  there  are  four  unknowns,  it  is  not  possible  to  directly  determine  the  spectrum 
parameters  from  ayy  and  &hh  alone,  so  further  assumptions  are  required  about  the 
relative  magnitudes  of  the  parameters  to  reduce  the  unknowns  to  two.  There  are 
two  methods  which  have  been  investigated. 


4.1  Method  1 


Beginning  with  data  for  the  horizontal  polarisation  cross  section  of  sea  clutter,  the  ratios 
of  the  Bragg,  whitecap  and  spike  terms  must  first  be  estimated,  and  then  they  can  be 
evaluated  from  ohh-  Assuming  equality  of  the  whitecap  component  in  both  polarisa¬ 
tions,  the  Bragg  component  of  the  vertical  polarisation  cross  section  can  be  calculated  by 
subtracting  the  whitecap  power  from  ayy. 

Treating  the  spike  component  Ps  in  Equation  55  first,  looking  upwind  at  a  wind  speed 
of  10  m/s,  Walker’s  data  indicates  that  the  spike  exceeds  the  horizontally  polarised  Bragg 
term  by  12.8  dB,  a  factor  of  19.1.  But  at  5  m/s,  spikes  are  rare  so  we  might  model  the 
ratio  to  diminish  to  0  at  4  m/s.  Also  we  note  that  sea  spikes  are  highly  directional  and 
do  not  exist  in  the  downwind  sector.  Thus,  they  might  be  modelled  as 


Ps 

PH 

rB 


=  Max 

[o,  19.1  LI 

<t> 

L  V  e  )\ 

COS  (j), 


=  0, 


0  <  (j)  <  7t/2, 
7t/2  <  <f>  <  IT 


(56) 


where  U  is  the  wind  strength  in  metres  per  second.  The  Max  function  is  used  to  set  the 
spike  component  to  0  for  wind  velocities  less  than  4  m/s,  and  the  cosine  function  reduces 
the  ratio  to  0  at  4>  =  7t/2.  To  reduce  the  ratio  more  quickly  with  aspect  we  could  consider 
using  the  second  or  a  higher  power  of  the  cosine  function. 

The  horizontally  polarised  Bragg  component  is  generally  small  compared  with  the 
whitecap  component  when  looking  in  the  upwind  direction,  but  downwind  these  compo¬ 
nents  are  roughly  equal.  Walker’s  upwind  data  indicates  that  at  10  m/s,  the  ratio  of  P\y 
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to  PjJ  is  11.5  dB,  or  a  factor  of  14.1.  Looking  downwind,  the  ratio  is  1.05  (0.2  dB)  at  7 
m/s  in  Walkers  data.  We  might  reasonably  assume,  in  the  absence  of  any  other  data,  that 
these  ratios  will  diminish  in  proportion  to  the  wind  velocity  since  wavebreaking  does  not 
occur  in  low  sea  states.  Finally  we  need  to  transition  the  ratio  from  its  upwind  value  to 
its  downwind  value  by  a  smooth  function,  and  thus  we  might  use 


Pw 

PH 

rB 


—  14.1  —  1 

(  1  +  COS  4>  \ 

|  +  1.05—  ( 

*  10  ' 

V  2  ) 

1  7  \ 

in  our  model. 


(57) 


Now  given  uhh  from  the  clutter  model,  we  can  determine  Pn  from  Equation  55,  thus 


PH  - 
DB  — 


&HH 


*  [Wb  +  Ww^  +  Ws^r 
B  B 


(58) 


The  values  of  Pw  and  Ps  then  immediately  follow  from  Equations  57  and  56.  Then 


rearranging  Equation  54,  Pg  can  be  determined  from 

y  _  °VV  —  y/nWwPw 


pv 

DB  — 


\TkWb 


(59) 


A  major  difficulty  with  this  model  is  that  near  broadside  the  power  in  the  whitecap 
component  \/n  WwPw  determined  from  (Thh  exceeds  ayy  at  higher  sea  states,  and  so  the 
predicted  value  of  Pg  is  negative.  The  solution  that  gives  acceptable  results  is  to  increase 
the  power  in  the  spike  component  near  broadside,  thus  reducing  the  whitecap  term.  This 
has  been  achieved  by  setting  the  spike  component  to  transition  to  0  at  (j)  =  n  instead  of 
7t/2,  thus 


Ps_ 

PH 

rB 


=  Max 


0,  19.1 


U  -  4 
6 


1  +  COS  (j) 


(60) 


in  the  model,  though  transitioning  to  0  at  0  =  3  7r/4  is  also  satisfactory.  It  might  be  argued 
that  spikes  do  not  exist  in  the  downwind  sector,  and  so  the  model  is  not  correct,  but  the 
problem  may  also  be  attributed  to  imprecision  in  determining  the  differences  between 
vertical  and  horizontal  polarisation  in  the  clutter  model(s).  Regardless,  the  total  power 
of  the  clutter  model  used  to  determine  the  spectral  components  is  preserved. 


4.2  Method  2 


In  this  method,  we  begin  the  analysis  with  ayv  to  determine  Pg  and  F\y  from  the  ratio 
of  these  terms,  and  then  determine  Pg  from  its  ratio  to  F\y ■  The  spike  term  is  then 
calculated  from  the  remainder  of  (Jhh  after  the  whitecap  and  horizontal  Bragg  components 
are  subtracted  from  it.  Walker  [26]  records  for  the  data  sets  he  analysed  that  upwind  at  a 
wind  speed  of  10  m/s,  the  ratio  of  F\y  to  Pg  is  0.437  (-3.6  dB),  and  downwind  at  7  m/s, 
the  ratio  is  0.0417  (-13.8  dB).  If  we  assume  that  these  ratios  are  proportional  to  the  wind 
speed,  i.e.  an  increased  wind  speed  increases  the  whitecap  component  proportionally,  and 
we  use  a  smooth  cosine  function  to  transition  from  upwind  to  downwind,  then 
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Pv 

rB 
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Rearranging  Equation  54 


pV  _ 


(Tvv 


WB  +  Ww^r 


(62) 


and  then  P\y  can  be  determined  from  Equation  61.  Likewise  PjJ  can  be  determined  from 
Equation  57;  it  is  interesting  to  note  that  no  information  about  (Jhh  is  used  in  this  method 
to  determine  the  horizontal  Bragg  component.  The  spike  term  then  follows  from 


Ps 


&hh  —  V^WbPb  ~  V^WwPw 
\PifWs 


(63) 


4.3  Scaling  the  Wavetank  Data  Model 

Open  ocean  spectral  data  appears  greater  in  frequency  and  wider  in  bandwidth  than  data 
obtained  from  wavetanks.  To  obtain  spectra  that  are  comparable  with  Walker’s  upwind 
and  downwind  examples,  the  models  for  the  frequencies  and  bandwidths  from  the  wavetank 
data  need  to  be  scaled  up.  In  doing  this  we  note  that  the  ocean  around  Portland  Bill  has 
strong  currents  which  may  influence  the  Doppler  frequencies,  particularly  the  frequency 
of  the  Bragg  components  of  vertical  and  horizontal  polarisation.  The  factors  that  were 
selected  by  visual  comparison  of  the  spectra,  and  are  used  to  modify  Equations  41  through 
52,  are 

1.  Bragg  peak  frequency  factor  -  1.3 

2.  Whitecap  peak  frequency  factor  -  2.7  looking  upwind,  2.3  looking  downwind 

3.  Bragg  bandwidth  factor  -  1.3 

4.  Whitecap  bandwidth  factor  -  2.7 

The  frequency  of  the  spike  component  is  the  same  as  the  whitecap  component,  and  its 
bandwidth  is  set  at  1.13  times  the  whitecap  bandwidth. 


5  Predictions  of  Sea  Clutter  Spectral 

Components 

Figures  9  to  12  show  each  of  the  clutter  spectral  components  obtained  from  the  three  mod¬ 
ified  clutter  models  discussed  earlier.  Both  parameter  extraction  methods  are  compared 
over  sea  states  1  to  5,  as  a  function  of  wind  aspect.  In  all  these  figures,  the  frequency  is 
9300  MHz,  and  the  grazing  angle  is  2°.  Note  that  for  the  spike  component  (Figure  12)  for 
method  1  only  sea  states  2  to  5  are  displayed  as  the  spike  component  is  set  to  0  for  wind 
strengths  less  than  4  m/s. 

The  horizontal  polarisation  Bragg  spectral  density,  (Figure  9),  shows  a  significant 
difference  between  the  two  methods.  In  method  2,  the  horizontal  polarisation  has  the 
bimodal  characteristics  of  ayv  although  ohh  is  unimodal  in  the  clutter  models.  The 
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Figure  9:  Bragg  component  of  horizontal 
(mf/mf/Hz),  expressed  in  dB,  as  a  function 
1  to  5,  at  a  grazing  angle  of  2° 
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Figure  10:  Bragg  component  of  vertical 
(rrif/mf/Hz),  expressed  in  dB,  as  a  function 
1  to  5,  at  a  grazing  angle  of  2° 
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Figure  11:  Whitecap  component  of  clutter  spectral  power  density  (ni2 /m? / Hz),  expressed 
in  dB,  as  a  function  of  wind  aspect  <f  at  9300  MHz  for  sea  states  1  to  5,  at  a  grazing  angle 

of  2° 
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Figure  12:  Spike  component  of  horizontal  polarisation  spectral  power  density 
(m2 /m2 / Hz),  expressed  in  dB,  as  a  function  of  wind  aspect  <f  at  9300  MHz  for  sea  states 
1  to  5,  at  a  grazing  angle  of  2° 
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explanation  for  this  is  that  in  this  method,  P[j  is  in  the  vicinity  of  15  dB  less  than  Pj?  for 
all  aspects  and  so  is  also  bimodal.  The  method  1  predictions  are  quite  different,  with  only 
about  a  2  dB  decrease  with  increasing  aspect  angle  until  the  downwind  sector  where  it  rises 
sharply  by  about  6  to  8  dB.  Even  for  sea  state  1,  for  which  there  is  no  spike  component, 
there  is  still  a  slight  peak  in  the  downwind  direction.  Even  though  the  horizontal  Bragg 
component  has  been  entirely  deduced  from  cthh,  which  is  unimodal,  the  peak  in  the 
downwind  direction  implies  a  bimodal  component  in  the  spectrum.  With  regard  to  the 
spacing  of  the  different  sea  states,  with  the  modified  GIT  model,  the  lines  for  sea  states 
2  to  5  are  tightly  bunched  together,  with  sea  state  1  quite  distinct  about  12  dB  less  in 
amplitude,  for  both  methods.  Note  also  that  with  the  modified  Hybrid  and  modified 
NRL  models  using  method  1,  the  lines  for  sea  states  1  and  2  almost  coincide  except  for 
downwind  aspects.  With  method  2,  the  spacing  of  the  Bragg  lines  obtained  from  the 
modified  hybrid  and  modified  NRL  models  corresponds  more  closely  to  the  spacing  in  the 
respective  clutter  models. 

The  vertical  polarisation  Bragg  component,  (Figure  10),  being  mostly  larger  than 
the  whitecap  component,  closely  follows  the  characteristics  of  ayy ■  For  all  three  clutter 
models,  the  variation  in  amplitude  with  aspect  for  method  1  is  greater  than  that  obtained 
with  method  2,  but  otherwise  quite  similar.  Bunching  of  the  curves  is  again  seen  for  the 
modified  GIT  model,  but  both  the  modified  Hybrid  and  modified  NRL  models  produce 
quite  regularly  spaced  sea  state  lines.  However  the  modified  NRL  model  with  method  1 
shows  unusual  behaviour  with  a  large  reduction  in  amplitude  for  higher  sea  states  near 
the  crosswind  look  direction,  resulting  in  the  sea  state  5  curve  crossing  both  sea  state  4 
and  3  curves,  and  almost  touching  the  sea  state  2  curve  in  this  region.  There  appears  to 
be  no  a  priori  reason  why  the  amplitude  of  the  Bragg  line  should  be  monotonic  with  sea 
state  in  this  region;  a  possible  explanation  could  be  that  the  larger  wave  action  of  sea  state 
5  with  associated  whitecaps  might  have  a  substantial  effect  on  propagation  of  the  Bragg 
component  in  the  crosswind  direction,  thus  resulting  in  the  dip  seen  with  method  1.  In 
the  absence  of  definitive  data  on  the  spectral  characteristics  of  the  clutter,  it  is  difficult 
to  decide  whether  the  predictions  on  method  1  with  the  modified  NRL  data  should  be 
rejected. 

The  Walker  model  of  sea  clutter  assumes  that  the  whitecap  component  of  sea  clutter 
is  identical  for  both  polarisations.  Beginning  first  with  the  development  of  the  whitecap 
component  from  cthh  of  method  1,  we  see  that  the  amplitude  is  quite  closely  unimodal 
except  of  a  slight  peak  in  the  higher  sea  states  in  the  downwind  direction.  The  predictions 
for  method  2  are  quite  different.  Initially  the  amplitude  decreases  at  the  rate  of  a  bimodal 
waveform  but  in  the  downwind  sector  is  almost  constant  from  about  100°  to  250°. 

In  method  1,  the  spike  component  (Figure  12)  is  forced  to  0  in  the  downwind  direction, 
but  in  method  2  is  allowed  to  account  for  the  remainder  of  (Jhh  after  determining  the  Bragg 
and  whitecap  components,  and  so  is  significantly  different;  in  the  latter  the  downwind 
component  is  about  12  dB  below  the  upwind  level  for  all  clutter  models.  It  is  questionable 
whether  a  spike  component  exists  in  the  downwind  sector,  but  it  is  necessary  to  include 
this  in  the  spectral  model  so  that  the  whitecap  component  is  the  same  for  both  vertical 
and  horizontal  polarisation  for  both  methods.  It  may  be  that  this  is  merely  a  problem 
inherent  in  the  averaging  of  the  clutter  data  in  each  polarisation  in  generating  the  clutter 
models  we  used  in  our  analysis. 
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When  we  look  at  the  total  spectra,  the  results  from  the  different  clutter  models  are 
quite  similar  in  form,  and  we  have  elected  to  show  only  the  results  obtained  with  the  mod¬ 
ified  NRL  clutter  model.  Figure  13  shows  the  combined  spectra  for  horizontal  polarisation 
for  both  methods  for  sea  states  1,  3  and  5  at  a  grazing  angle  of  2°  for  wind  aspects  in  45° 
increments  between  upwind  and  downwind.  Reasonable  correspondence  is  obtained  with 
Walker’s  measurements.  The  Bragg  term  is  relatively  insignificant  and  mainly  manifests 
itself  as  a  slight  bump  in  the  lower  half  of  each  spectrum  curve.  In  the  vertical  polarisa¬ 
tion  curves  (Figure  14)  the  Bragg  term  is  dominant,  especially  at  low  sea  states.  Close 
correspondence  with  Walker’s  measurements  is  achieved  with  a  sea  state  of  4.22,  for  which 
the  wind  velocity  is  10  m/s,  according  to  the  relationship  between  wind  velocity  and  sea 
state  used  in  the  clutter  model. 

The  differences  between  the  spectra  obtained  with  the  two  methods  are  relatively  minor 
but  it  is  not  known  whether  they  are  significant  in  detection  calculations. 


6  Conclusions 

This  investigation  has  demonstrated  that  it  is  possible  to  deduce  a  credible  spectrum  for 
sea  clutter  from  standard  total  cross  section  sea  clutter  models,  with  the  support  of  a 
set  of  assumptions  about  the  relative  relationship  of  the  components  that  comprise  the 
sea  clutter  spectrum.  However,  it  has  generated  as  many  questions  as  it  has  answered; 
all  the  results  appear  acceptable,  and  the  different  generating  clutter  models  will  need  to 
be  used  advisedly  for  the  purpose  to  which  the  spectrum  data  is  to  be  put.  In  carrying 
out  this  investigation,  the  modified  NRL  model  has  some  features  which  might  make  it 
more  suitable  for  this  use;  in  particular,  the  vertical  and  horizontal  polarisation  terms 
are  separately  optimised  to  conform  with  Nathanson’s  data,  not  relying  upon  a  fixed 
dB  difference  independent  of  sea  state,  as  used  in  the  other  two  models,  which  may  be 
significant  in  determining  the  spectral  components.  Also  the  much  improved  accuracy  at 
low  sea  states  is  significant  for  the  requirement  for  which  the  model  was  developed. 

An  experimental  measurement  program  which  investigates  the  validity  of  the  assump¬ 
tions  in  deriving  these  results  is  warranted  as  there  is  a  need  for  accurate  low  grazing  angle 
Doppler  sea  clutter  models  across  a  number  of  Defence  programs.  In  particular  there  is 
need  for  the  measurement  of  the  separate  terms  comprising  the  spectrum,  to  determine 
their  ratio  at  various  wind  strengths  and  a  wide  range  of  aspect  angles,  as  the  predictions 
of  this  research  depend  on  these  ratios.  It  is  also  desirable  to  confirm  Walker’s  model, 
that  the  whitecap  component  is  identical  for  both  polarisations,  for  all  aspects,  as  it  does 
not  seem  justified  from  an  electromagnetic  scattering  perspective.  Data  is  also  required 
on  the  behaviour  of  the  spike  component,  particularly  its  dependence  on  the  wind  aspect 
angle  and  whether  it  is  present  in  the  downwind  sector. 

It  is  hoped  that  this  model  will  prove  useful  for  more  applications  than  the  specific 
project  which  it  was  developed,  as  it  is  an  appropriate  tool  to  assist  in  the  assessment 
of  performance  of  a  number  of  radar  systems  of  current  interest  to  Defence.  Given  the 
highly  speculative  nature  of  the  development  of  this  model,  it  is  anticipated  that  the 
assumptions  and  the  methods  for  deducing  the  spectral  components  will  evolve  as  new 
data  comes  available. 
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Figure  13:  Predicted  power  spectral  density  (m2 /m2  /  Hz),  expressed  in  dB,  of  horizontal 
polarisation  sea  clutter  at  9300  MHz,  using  the  modified  NRL  clutter  model,  at  a  grazing 
angle  of  2° ,  for  wind  aspect  f)  between  0°  and  180° 
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Figure  1 f:  Predicted  power  spectral  density  (rnf/mf/Hz),  expressed  in  dB,  of  vertical 
polarisation  sea  clutter  at  9300  MHz,  using  the  modified  NRL  clutter  model,  at  a  grazing 
angle  of  2° ,  for  wind  aspect  $  between  0°  and  180° 
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